Structure of Thermal Quasifermion in the QCD/QED Medium 
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In this paper we carried out a nonperturbative analysis of a thermal quasifermion in the chiral 
symmetric thermal QCD/QED medium by studying its self-energy function through the Dyson- 
Schwinger equation (DSE) with the HTL resummed improved ladder kernel. 

Our analysis reveals several interesting results, some of which may force us to change the image 
of thermal quasifermions; (1) The thermal mass of a quasifermion begins to decrease as the strength 
of the coupling gets stronger and finally disappears in the strong coupling region, thus showing a 
property of a massless particle, and (2) its imaginary part (i.e., the decay width) persists to have 
0(g 2 Tlog(l/g)) behavior. These results suggest that in the recently produced sQGP, the thermal 
mass of a quasiquark should vanish. Taking into account the largeness of the imaginary part, it 
seems very hard for a quark to exist as a qausiparticle in the sQGP phase. 

Other important findings are as follows; (3) The collective plasmino mode disappears also in the 
strongly coupled system, and (4) there exisits an ultrasoft third peak in the quasifermion spectral 
density at least in the weakly coupled QED/QCD plasma, indicating the existence of the ultrasoft 
fermionic mode. 



I. INTRODUCTION 

It is believed that the Relativistic Heavy Ion Collider 
(RHIC) at BNL and the Large Hadron Collider (LHC) 
at CERN have produced the primordial state of matter, 
namely the quark-gluon plasma (QGP), and liberated the 
quark and gluon degrees of freedom. Subsequent analyses 
have shown that the produced QGP medium shows the 
property close to that of a perfect fluid. This fact leads 
us to the understanding that the QGP produced in the 
energy region of the RHIC and LHC is a strongly inter- 
acting system of quarks and gluons, namely, the strongly 
coupled QGP (sQGP) [j. ' 

Since the discovery of sQGP phase, the behavior and 
the properties of quasiquark in the new sQGP phase have 
attracted much attention; does the quark still work as the 
basic degree of freedom in the new phase or not ? It is 
also pointed out theoretically that the hadronic excita- 
tion affects the spectral density of quasiquark even in the 
chiral symmetric phase, thus showing some characteristic 
structures near the phase boundary |2|. 

Up to now, most of the theoretical findings on thermal 
quasiquarks in the QGP are obtained through analyses 
with the assumption of weakly coupled QGP at high tem- 
perature, i.e., analyses through the hard-thermal- loop 
(HTL) resummed effective perturbation calculation Q, 
or those through the one-loop calculation with the mas- 
sive bosonic mode, or by replacing the thermal gluon with 
the massive vector boson Kitazawa et al. Q have 
pointed out the three-peak structure of the quasifermion 
spectral density and the existence of the massless third 
mode. Such analyses, however, cannot be justified in 
studying the thermal quasiparticle in the sQGP created 
in the energy region of RHIC. What we need is the non- 
perturbative analysis to explore the the properties of 
strongly coupled system. 

Nonperturbative calculations of correlators within lat- 
tice QCD are performed in Euclidean space and give in- 



teresting results [5]. However, strictly speaking it is not 
possible to carry out an analytic continuation that is nec- 
essary to determine the spectral function. In addition, it 
is difficult on the lattice to respect the chiral symmetry 
that should be restored in the sQGP phase, though we 
are interested in the property of thermal quasiparticlcs 
in the chiral symmetric sQGP phase. 

In this paper we perform a nonperturbative analysis of 
a thermal quasifermion in thermal QCD/QED by study- 
ing its self-energy function through the Dyson- Schwinger 
equation (DSE) with the HTL resummed improved lad- 
der kernel. Our analysis may overcome the problems in 
the previous analyses listed above for the following rea- 
sons: (1) it is a nonperturbative QCD/QED analysis, (2) 
we study the DSE in the real-time formalism of thermal 
field theory, which is suitable for the direct calculation 
of the propagator, or the spectral function, (3) we use 
the HTL resummed thermal gauge boson (gluon/photon) 
propagator as an interaction kernel of the DSE, and take 
into account the quasiparticle decay processes by ac- 
curately studying the imaginary part of the self-energy 
function, and finally, (4) we present an analysis based on 
the DSE that respects the chiral symmetry and describes 
its dynamical breaking and restoration. Our analysis is 
nothing but an application of our formalism employing 
the DSE to the study of thermal quasifermions on the 
strongly coupled QCD/QED medium with chiral sym- 
metry 0, 0]- 

With the solution of the DSE with the HTL resummed 
improved ladder kernel, we study the properties of the 
thermal quasifermion spectral density and its peak struc- 
ture, the dispersion law of the physical modes corre- 
sponding to the poles of thermal quasifermion propaga- 
tor, through which we elucidate the properties of thermal 
mass and the decay width of fermion and plasmino modes 
and also pay attention to properties of the possible third 
mode, both especially in the sQGP phase. 

Analogous studies employing the DSE are carried out 
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by several groups Q. All these analyses solve the DSE 
in the imaginary-time formalism, and try to perform an 
analytic continuation. Harada et al. study the DSE with 
a ladder kernel in which the tree-level gauge boson prop- 
agator is used, while Qin et al. and Mueller et al. use the 
maximum entropy method to compute the quark spectral 
density. Qin et al. also paid a special attention to the 
massless third mode. 

Our analysis reveals several interesting results, some 
of which may force us to change the image of thermal 
quasifermions: (1) The thermal mass of a quasifermion 
begins to decrease as the strength of the coupling gets 
stronger and finally disappears in the strong coupling 
region, thus showing a property of a massless particle, 
and (2) its imaginary part (i.e., the decay width) persists 
to have 0(g 2 T\og(l/g)) behavior. These results suggest 
that in the recently produced sQGP, the thermal mass 
of a quasiquark should vanish. Taking into account the 
largeness of the imaginary part (i.e., the decay width), it 
seems very hard for a quark to exist as a qausiparticle in 
the sQGP phase. 

Focusing on the fact (1) above, we have already re- 
ported briefly in Rcf. [9], and in the present paper we 
give results of detailed analysis. The fact (4) has also 
been pointed out briefly in Ref. [9J, and will be studied 
further in detail in a separate paper [Io| . 

Other important findings are as follows; (3) The collec- 
tive plasmino mode disappears also in the strongly cou- 
pled system, and (4) there exisits an ultrasoft third peak 
in the quasifermion spectral density at least in the weakly 
coupled QED/QCD plasma, indicating the existence of 
the ultrasoft fermionic mode. 

This paper is organized as follows; In Sec. HJ we 
present the HTL resummed improved ladder Dyson- 
Schwinger equation for the quasifermion self-energy func- 
tion, with which we investigate the property of ther- 
mal quasifermion in the chirally symmetric QGP phase, 
and give the results in Sec. IIIII In Sec. MI A| properties 
of the quasifermion spectral density is studied, and in 
Sec. MI Bl we discuss the problem in the relation between 
the peak-position of spectral density and the zero-point 
of the inverse quasifermion propagator. The dispersion 
law of the quasifermion is studied in Sec. MI CI and the 
vanishment of thermal mass and the disappearance of 
the plasmino mode in the strongly coupled system are 
pointed out. Properties of the thermal mass and the 
existence of the third peak or the ultrasoft mode are dis- 
cussed in Sec. MI Dl and Sec. MI El respectively. Finally in 
Sec. IIIIF1 properties of the decay width of quasifermion 
is studied. Summary of the paper and discussion is given 
in Sec. IIVI Several appendices are also given. In Ap- 
pendix |A] we explain the approximations to get the HTL 
resumed improved ladder DSE to be solved. The cutoff 
dependence of our analysis is discussed in Appendix [Bj 
and the phase boundary between the chirally symmetric 
and broken phases in the Landau gauge is briefly ex- 
plained in Appendix [C] Finally Appendix [D] is devoted 
to explain why we do not use the peak position of the 



spectral density as the condition to determine the parti- 
cle onshellness. 



II. THE HTL RESUMMED IMPROVED 
LADDER DYSON-SCHWINGER EQUATION 

In this paper we study the thermal QCD/QED in 
the real time closed time-path formalism 11], and solve 
the Dyson-Schwinger equation (DSE) for the retarded 
fermion self-energy function E^j, to investigate the prop- 
erty of thermal quasifermion in the chiral symmetric 
QGP phase. Throughout this paper we study the mass- 
less QCD/QED in the Landau gauge. 

As is well-known, at zero-temperature the Landau 
gauge plays an essential role to ensure the gauge invari- 
ance of the solution of the ladder DSE because it is proved 
that in the Landau gauge the wave-function receives no 
renormalization, i.e., A(P) = 1 [l2l [j^. At finite tem- 
perature, however, A(P) ^ 1 even in the Landau gauge, 
and there is no special reason to choose Landau gauge 
any more. In this analysis we choose the Landau gauge 
firstly for the sake of simplicity, and secondly for the sake 
of comparison with other works. 

In this section we present the HTL resummed DSE for 
and also give an explication about the improved lad- 
der approximation we make use of to the HTL resummed 
gauge boson propagator. We also calculate the effective 
potential for the retarded fermion propagator Sr in or- 
der to find the "true solution" when we get several "so- 
lutions" of the DSE. 



A. HTL resummed improved ladder DSE for 
fermion self-energy function Eij 

The retarded quasifermion propagator Sr(P),P 
(po,p), is expressed by 



Sr(P) 



P+ ie 7 ° - Efl(P) ■ 



(1) 



The retarded fermion self-energy function E^ can be 
tensor-decomposed in a chiral symmetric phase at finite 
temperature as follows; 



Efl(P) = (1 - A(P)) Pl Y - P(P)7 C 



(2) 



A(P) is the inverse of the fermion wave function renor- 
malization function, and B(P) is the chiral invariant 
mass function. The c-number mass function does not 
appear in the chiral symmetric phase. 

In the real time closed time-path formalism, by adopt- 
ing the tree vertex and the HTL resummed gauge bo- 
son propagator for the interaction kernel of the DSE, 
we obtain, in the massless thermal QED/QCD, the HTL 
resummed improved ladder DSE for retarded fermion 
self-energy function E# [(| [l4| (coupling a = g 2 /iTT : 
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g 2 = g 2 s C f for QCD, g = e for QED) 
g 2 f d A K 

x [*r^(-P,ir,p-^)^(-x,x) 

x*r^(-^, P, A - PYGbj^^K -P,P-K) 
+* t raa(- p > K ' p - K)S rr (-K, K) 
x*r U AAR(-K, P,K- P)*G RA ^{K -P,P-K)]. 

(3) 

Here * G, lu is the HTL resummed gauge boson propaga- 
tor [T5], [l|| where P = RA and C = PP denote the re- 
tarded and the correlation components, respectively, and 
= 7^ in the present approximation. 
There are many attempt to carry out the higher order 
calculation within the HTL resummed effective pertur- 



bation theory and to get information beyond the appli- 
cability region of the HTL approximation 0, [l7], [Hj] . The 
DSE with the HTL resummed gauge boson propagator 
as an interaction kernel can take the dominant effects 
of thermal fluctuation of O(gT) into account nonpertur- 
batively. Thus we expect the HTL resummed improved 
ladder DSE to enable us to study wider regions of the 
couplings and temperatures, e.g., the strongly coupled 
QGP medium, than those restricted by the HTL approx- 
imation, i.e., the regions of weak couplings and high tem- 
peratures. 

The explicit expression of the HTL resummed im- 
proved ladder DSE in the Landau gauge to determine 
the scalar invariants A(P) and B(P) in the chiral sym- 
metric QGP phase becomes coupled integral equations as 
follows; 



p 2 [l-A(P)] 



d K 



(2tt)4 
+2p k g< T og P o} 



{1 + 2n B (po - fco)}Im[ *G%(P - K)\ x 



{K a P p + K p P a - p {K a g p0 + K p g a0 ) - k (P a g p0 + P p g a0 ) + pkzg ap 

A(K) 



-2poS<xOffpo} 



[k + B{K) + ie} 2 - A{K) 2 k 2 
k + B{K) 



[feo + B(K) + ie} 2 — A(K) 2 k 2 



{P<r9p0 + Pp9crO 

+ {1 - 2n F (k )} x 



f G^(P - K)Im {K a P p + K p P a - p (K a g p0 + K p g a0 ) - k (P a g pQ + P p9<j0 ) 



+pkzg ap + 2p k gcrog p o} 



A{K) 



{Pa9p0 + P p g<rO - 2pog<rOg P o} 

d 4 K 



[ko + B(K) + ie} 2 - A(K) 2 k 2 
h + B(K) 



[k + B(K)+ie] 2 - A{K) 2 k 2 
{1 + 2n B ( P o - fco)}Im[ *G P °(P - K)\ x 

A(K) 



(2tt)* 

{Kagpo + K p g aQ - 2k g a0 g p0 } 



-{2g p0 2g aQ - g ap \ 
*G P °(P-K)lm 
-2k Q g a0 g p0 } 



[k a + B(K) + ie} 2 - A{K) 2 k 2 
fco + B{K) 



[feo + B(K) + ie} 2 - A(K) 2 k 2 
A{K) 



[k + B(K)+te] 2 -A(K) 2 k 2 
k + B{K) 



{l-2n F (fc )} x 
{ K a g p o + K p g a0 



{2g p0 2g a0 ~ gap} 



(4a) 



(4b) 



[k + B(K)+ie] 2 - A(K) 2 k 2 
where n B {x) and uf(x) are the Bosc-Einstein and the Fermi-Dirac equilibrium distribution functions, respectively, 

1 1 



n B (x) 



exp(a;/T) - 1 



exp(x/T) + 1 



r 



The above DSEs, Eq. (jj), are still very tough to be We thus adopt the instantaneous exchange approxima- 
attacked, and need further approximations to be solved. tion to the longitudinal gauge boson propagator, i.e., we 
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set the zero-th component of the longitudinal gauge bo- 
son momentum q to zero. Details of the approximation 
we use are explained in Appendix [A] 

In solving the DSEs, Eq. we are forced to intro- 
duce a momentum cutoff in the integration over the four- 
momentum j d A K, K = (&o,k). We use the following 
cutoff method (A denotes an arbitrary cutoff parameter 
and plays a role to scale any dimensionful quantity, e.g., 
T = 0.3 means T = 0.3A); 

three momentum k : A; — |k| < A 

energy k : \k \ < A (A = A - 5A) 

We determine the parameter Ao so as to get a stable so- 
lution for the fermion spectral density. Over the range of 



the temperature and the coupling we study in the present 
analysis, we set A = 2A. In fact solution is totally stable 
for A > 2A. In Appendix IBI we give details of the cutoff 
dependence of the solution. 



B. The effective potential V[<Sr] for the retarded 
full fermion propagator Sr 

Above DSEs, Eq. (T5|), may have several solutions, and 
we choose the "true" solution by evaluating the effective 
potential V^Sj?] for the fermion propagator function S R , 
then finding the lowest energy solution. The effective 
potential is expressed as [19( 



V{S R }=Tr{fS R }+Tr\n[S B 1 }+^ J^J [l^ R {K) lv S R {P)G^ {P - K) 

+ 1 »Sc(K) lv S R {P)G R u {P - K) + lfl S R (K) lu Sc(P)G%(P - K)] , (5) 

I 



where the 1st and the 2nd terms correspond to the one- 
loop effective potential, while the 3rd term corresponds 
to the two-loop contribution. 



III. SOLUTION OF THE HTL RESUMMED 
IMPROVED LADDER DYSON-SCHWINGER 
EQUATION 

In this section we give the solution of the HTL re- 
summed improved ladder DSE for the retarded self- 
energy function S_r, and study its consequences in the 
chiral symmetric phase. Special attention is paid to 
the consequences in the strongly coupled QCD/QED 
medium, in order to get an information on the ther- 
mal quasifermion in the strongly coupled QGP (sQGP) 
phase, recently discovered through the experiments at 
RHIC and LHC [1]. Part of the results are already briefly 
reported @- 



A. quasifermion spectral density 

1. Spectral density of quasifermion p± 

In the chiral symmetric QCD/QED phase, the 
quasifermion propagator can be expressed as 



with 



S R (P) 
where 



1 

D~4 



7 o + P_a 
p 



i 

d7 



o _ PO_ 
P 



D ± (P)=p + B(p ,p)TpA{p ,p) 



(6) 



(7) 



Re[D + (p 0)P )} = -Re[D_(-p 0>P )}, (8a) 
lm[D + (p ,p)] = lm[D_(-p ,p)}. (8b) 

The spectral density of quasifermion p± is defined by 
P±(po,p) = --Im- 



7T D±(P) 

-Ilm- 



1 



7T p + B(p ,p) tpMpo,p) : 
which satisfies the sum rules [2(| 

dpoP±(Po,p) = 1, 

/OO 
dpoPoP±(po,p) = ±P, 
-CO 

/OO 
dpoPoP±(po,p) = P 
-OO 

It also satisfies the symmetry property 
P±(po,p) = Pt(~Po,p) 



2 i 2 

1 mi. 



(9) 

(10a) 
(10b) 
(10c) 

(11) 



We therefore study only p + (po,p) throughout this paper. 

To study the spectral properties of quasifermion in the 
chiral symmetric QGP phase, we must at first make sure 
that we are in fact studying inside the chiral symmetric 
phase. We have already studied the phase structure of 
thermal QCD/QED through the same DSE approach 0, 
and determined the phase boundary between the chiral 
symmetric and the chiral symmetry broken phases. Thus 
we are sure that the region of temperatures and couplings 
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we study in the present paper are well within the chiral 
symmetric phase. For details, see Appendix [CI where we 
give the results of our analyses to determine the phase 
boundary. 

It is to be noted that the temperature T in the present 
analysis represents the temperature scaled by the mo- 
mentum cutoff A, thus < T < 1, and it does not denote 
the real temperature itself. With the temperature thus 
defined, we can determine the phase boundary curve in 
the two-dimensional (a, T)-plane, separating the chiral 
symmetry broken/restored phase [f|. If we measure the 
temperature T relative to the critical temperature T c , 
high or low temperature has a definite meaning 



2. Structure of quasifermion spectral density 

Now let us present the properties of quasifermion 
spectral density p + computed by the solution of DSEs, 
Eq. ([4]). At first we should note the fact that the 
quasifermion spectral density p + thus determined well 
satisfies the sum rules Eqs. (|10a|) and (|10b| within a few 
present error. This fact proves a posteriori the adequacy 
of our choice of cutoff parameter Ao = 2A in the region 
of the couplings and the temperatures we study. The 
third sum rule is heavily depends on the HTL calcula- 
tion, and the agreement depends on the couplings and 
the temperatures. 

Then let us show the structure of p+ (po , p) &s a func- 
tion of po and p in the two-dimensional (po,p)-pl&ne. 
For convenience we study the dimensionless quantity 
P-\-(jPo,p)m*p where m*j denotes the thermal mass deter- 
mined through the next-to-leading order calculation of 
the HTL resummed effective perturbation theory [H, [2l[ , 



mf 



i- A JL 



9_ 

'2tt 



9_ 

4tt 2 



(12) 



In measuring at moderately high temperature T = 0.4, 
we can see in Fig. Q] the three typical peak-structures 
depending on the strength of the coupling a = g 2 /Ait] 

i) At weak coupling a = 0.005 we can see three peaks as 

a function of po at p — 0. Two sharp peaks of them 
at positive and negative po represent the fermion 
and the collective plasmino modes, respectively (T5L 
HH , and the slight third "peak" barely recognizable 
around po — corresponds to the massless, or the 
ultrasoft mode [1, [22| . The plasmino mode and the 
massless mode rapidly decrease and disappear as 
the size of momentum p becomes large. 

ii) At the intermediate strength a — 0.2, we can see 
only two peaks at pg ^ as a function of po even 
at p = 0, and unable to recognize the existence of 
the third peak corresponding to the massless pole 



in this region of the coupling. The peak at the 
negative side of po axis that may corresponds to 
the collective plasmino pole rapidly disappears as 
p gets large. 

iii) At the strong coupling a — 1.0 we can only recog- 
nize, at any size of the momentum p, the existence 
of a broad "peak" that may represent the massless 
pole. No massive pole exists in the strong coupling 
region. 

It should be noted the vast differences of the height 
of the peak and the spread of the spectral density in 
the three cases of the coupling strength i), ii) and iii), 
clearly showing the broadness of the "peak" in the strong 
coupling environment. 

Figure [T] shows the coupling a-dependence of the 
quasifermion spectral density p+m*j at fixed temperature 
T = 0.3. We can also see the temperature T-dependence, 
which is given in Fig. [5] 

In Fig. [TJ as noted above, we see the transition of the 
peak structure of spectral density as the strength of the 
coupling varies with the temperature kept fixed: triple 
peaks at small couplings, double peaks at intermediate 
couplings and finally single peak at strong couplings. 

Figure [2] shows that the analogous behavior is also ob- 
served when the temperature of the environment varies 
with the strength of the coupling kept fixed. At small 
couplings (a = 0.01), the three-peak structure at low 
temperature (T = 0.2) changes to the double-peak struc- 
ture at high temperature (T = 0.8), and at intermediate 
couplings (a — 0.45), the double-peak structure at low 
temperature (T = 0.2) tends to the single-peak structure 
at high temperature (T = 0.8). 

Here let us see more carefully the structure of spectral 
density at p = 0, p+(po,p = 0), as a function of p$. The 
Po-coordinate of the peak position of p + (po,p = 0) will 
give the mass of the corresponding mode. In Fig.|3]shown 
are the spectral densities, p+(poiP = 0), at moderately 
high temperature T = 0.3, in the weak coupling region 
a < 0.01, in the region of intermediate coupling strength 
a « 0.1 ~ 0.2, and in the strong coupling region a ~ 1. 

In order to see what actually happens during the tran- 
sition from the triple peak structure in the weak coupling 
region to the double peak one in the region of interme- 
diate coupling strength, and finally to the single peak 
one at strong couplings, we present in Figs. @] and [5] the 
Re[D+(p ,p = 0)] and lm[B(p ,p = 0)] at T = 0.3, re- 
spectively, both of which show three curves correspond- 
ing to the three regions of the coupling as in Fig. [3] 

At weak coupling a — 0.005, we can clearly see in Fig. [3] 
two sharp peaks at positive and negative po, representing 
the quasifermion and the plasmino poles. The center- 
positions of both of peaks are in fact at po/m*^ ~ ±1, 
which in fact almost coincide with the solution of the on- 
shell condition Re[D + (p ,p = 0)] = 0, as can be easily 
seen in Fig. 01 Thus at the weak coupling and high tem- 
perature both the quasifermion and the plasmino modes 
have a common thermal mass mZ, Eq. I12[ which is, as 
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FIG. 1: The coupling a-dependence of the quasifermion spectral density p+(po,p) at T = 0.3. 




FIG. 2: The temperature T-dependence of the quasifermion spectral density p+(po,p) at the weak coupling a = 0.01 (upper 2 
graphs) and at the strong coupling a = 0.45 (lower 2 graphs). 



already noted, determined through the next-to-leading 
order calculation of HTL resummed effective perturba- 
tion theory @,[2l|. 

We can also barely recognize the existence of a slight 
"peak" around p = 0, corresponding to a massless pole. 
The existence of this massless mode is also indicated by 
the onshell condition Re[D + (po,p = 0)] = 0, in which 
Pa = is always a solution. The problem on this third 
"peak" will be discussed in the Sec. MI El 

It is also worth noticing that at weak coupling a = 
0.005 Im[B(po,p = 0)] shows a three peak structure; a 
sharp steep peak centered at p — 0, and two slight peaks 
centered at \po\ = m*^ (23J. All three peaks in Im[B] ap- 
pear corresponding to the onshell point Ke[D + (po,p = 



0)] = 0, as were the case in the peaks of the spectral den- 
sity, the sharpness of the peaks being completely turned 
over. 

At intermediate coupling a = 0.2, the fermion spectral 
density exhibits a typical double peak structure. The 
existence of these peaks, however, is not easy to under- 
stand. As we can easily make sure by comparing Fig. [3] 
and Fig. [H they do not have exact correspondence to the 
poles of the fermion propagator, i.e., the zero-point of the 
inverse propagator Rc[D + (p ,p = 0)] = 0. 

At strong coupling a — 1.0, the situation becomes 
very simple. There exists only a broad single peak, 
whose existence is indicated by the onshell condition 
Ke[D + (po,p = 0)] = 0, in which po — is the only 
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1.4 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 
P /m *f 

FIG. 3: The spectral density at p — 0, a+(po,P = 0), at T = 
0.3, in the weak coupling region a = 0.05, in the intermediate 
coupling region a — 0.20 and in the strong coupling region 
a = 1.0. 
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FIG. 4: The real part of the inverse fermion propagator at 
p = 0, Re[D+(po,P = 0)], at T = 0.3, in the weak coupling 
region a — 0.05, in the intermediate coupling region a — 0.20 
and in the strong coupling region a — 1.0. 
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FIG. 5: The imaginary part of the inverse fermion propagator 
at p = 0, lm[D+(p ,p = 0)]=lm[B(p o ,p = 0)], at T = 0.3, 
in the weak coupling region a = 0.05, in the intermediate 
coupling region a = 0.20 and in the strong coupling region 
a = 1.0. 



solution in the strong coupling region, see Fig. [4] It is 
not clear whether the massless peak at strong couplings 
is exactly the same one at weak couplings noted above, or 
not. We will discuss on this massless mode in Sec. IIIIEI 

To understand the typical structure of the fermion 
spectral density explained above, it is always important 
to correctly take notice of the height of the peak and 
the width of the corresponding pole, namely the height 
of the peak and the width of lm[B(P)}, given in Fig. [SJ 
in connection with the solution of the onshell condition 
Re[£)+(p o ,P = 0)] =0. 

In this sense the appearance of the double-peak struc- 
ture at intermediate couplings is somewhat confusing. It 
is because, while there is a clear correspondence between 
the triple peaks at small couplings and the physical poles 
or modes (i.e., the quasifermion, the plasmino and the 
massless or ultrasoft modes), the double peaks at inter- 
mediate couplings do not have an obvious correspondence 
to the physical poles or modes. They may correspond 
to the quasi-fermion and the plasmino modes, but the 
center-positions of the peaks are apparently bigger than 
expected from the value of thermal mass, see, Fig. [3] 
In addition, as explained above, the third peak corre- 
sponding to the massless or the ultrasoft mode can not 
be recognized at all. This problem might arise from the 
broad-peak structure of the imaginary part of the mass 
function, Im[B(p ,p)], centered at po = 0, see, Fig. [SJ 
and will be discussed later in Sec IIIIB1 Sec. MI CI and 
Appendix |DJ 

With the appearance of this problem, it seems better 
to determine the position of the quasi-fermion pole by 
the solution of the onshell condition Re[D + (po,p)} = 
than to determine by the peak-position of the spectral 
density. 



B. What determines the peak-position of spectral 
density? Or, the relation between the peak-position 
of spectral density p+(P) and the zero-point of the 
inverse propagator Re[D+(P)] 

Here we study the problem; What determines the 
peak-position of spectral density? At the end of the last 
Sec lIII A~2l we have briefly commented on the problem 
by focusing on the relation between the peak-position 
of spectral density p+(P) and the zero-point of the in- 
verse propagator Re [D + (P)]. There we also noticed that 
we should correctly take into account the information on 
the Im[P(P)]. 

Let us summarize what we have disclosed, a) At 
weak couplings there are two sharp peaks located at 
Po/mJ: ~ ±1, which in fact almost coincide with the so- 
lution of the onshell condition Re[D + (p ,p = 0)] = 0. 
The third slight "peak" around po = is also indicated 
by the onshell condition Re[D + (p ,p = 0)] = 0, in which 
Po = is always a solution, b) Typical double peaks at 
intermediate strength of coupling do not have an obvious 
correspondence to the zero-point of the inverse propaga- 
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tor Ke[D + (p ,p = 0)] = 0. c) At strong couplings, there 
exists only one broad "peak", whose existence is indi- 
cated by the onshell condition Re[D + (p ,p = 0)] = 0, in 
which po = is the only solution in the strong coupling 
region. 

There appear two questions. 1) What happens at in- 
termediate strength of coupling? 2) What causes the 
huge difference of the peak height between sharp peaks 
of fermion and plasmino modes and the slight peak of 
massless mode at weak couplings? 

We can add one more question; Does the massless peak 
(or the pole) at strong couplings represent the same mass- 
less mode at weak couplings? This third question, how- 
ever, will be discussed in a separate paper. 

Now let us study questions 1) and 2) in order. 

On the question 1): First let us see the solution of the 
onshell condition Re[D + (p ,p = 0)] = 0. In the range of 
intermediate couplings around a ~ 0.2 (temperature is 
fixed at T — 0.3), the real part of the chiral symmetric 
mass function at p = 0, Re[B{po,p = 0)], as a function 
of po exhibits a subtle structure around the origin. In 
studying the small po region, it has a steep valley /peak 
structure at weak couplings, but as the coupling becomes 
stronger this valley /peak structure eventually diminishes 
in size and begins to behave almost as a straight line. 

Intermediate coupling region is the transition region; 
As the coupling gets stronger the two solutions of the 
onshell condition at |po| 7^ eventually approach to 
Po = and coincide with the solution at po — that 
always exists irrespective of the strength of the coupling. 
Thus the number of solutions of the onshell condition 
Re[D + (po,p — 0)] =0 changes suddenly from three to 
one. 

We should check here, in the considered region at 
T = 0.3 with couplings around a = 0.2 and stronger, 
where the real part of the inverse propagator vanishes, 
i.e., where the solutions of Re[D + (po,p = 0)] = exist. 
There are three solutions, two of them sit at po ^ 0, i.e., 
\p \ — 0.6 and the third one at p = 0, see Fig. [U thus 
indicating the existence of three poles, or the appearance 
of three peaks in the spectral density. 

We should then see the shape and the position of the 
peak of lm[B(po,p = 0)]. lm[B(po,p = 0)] always has 
a single broad peak around po — in the corresponding 
region of temperatures and couplings, i.e., T — 0.3 and 
couplings around a — 0.2 and stronger, see Fig. [5j 

With these fact we understand that, at intermediate 
coupling a = 0.2, the peak structure of lm[B(p ,p = 0)] 
plays an important role, scratching out (washing away) 
the peak of the spectral density at po — 0, and the 
not-so-steep but still gaussian decreasing structure of 
lm[B(po,p — 0)] makes the positions of the peaks of the 
spectral density at \po\ ~ 0.6 shift to larger \po\ values. 

On the question 2): To understand this question, let 
us see Figs. [3l HI and [5] in the weak coupling region. The 
spectral density exhibits two sharp peaks at po ~ ^"^/j 
and one barely slight peak at po — 0. The positions 
of these three peaks exactly agree with the three solu- 



tions of the onshell condition Re[D + (po,p = 0)] = at 
Po ~ and at po = 0, thus these three mode rigidly 

correspond to the fermion, plasmino and massless modes, 
respectively 0, GJ, EE SI . 

In contrast with the spectral density, the structure 
of Im[B(po,p = 0)] (—Im[D + (p ,p — 0)]) is simple. 
lm[B(po,p — 0)] at weak coupling a = 0.005, as can 
be seen in Fig. [5j exhibits a sharp peak at po = and 
two slight peaks at |po|/w1 ~ 1. These peaks have clear 
correspondence to the three solutions of the onshell con- 
dition Re[D + (p 0} p = 0)] = 0. At the positions of two 
sharp peaks in the spectral density it is essential that 
Re[Z? + ] is zero, and the imaginary part of it, or Im[£>], 
is so small that it does not play any essential role in 
the structure of the spectral density. At the position of 
massless pole po = 0, however, Im[B] is so large that it 
plays an important role to almost scratch out the fact 
that Re[Z? + ] is zero, thus the peak structure almost dis- 
appears at around the origin. 

The peak height at the pole is determined by its pole- 
residue. This fact means that, by measuring the ratio 
of peak-heights between the sharp peak representing the 
quasifermion mode and the slight peak at the origin rep- 
resenting the massless or ultrasoft mode, we can deter- 
mine the ratio between the corresponding pole-residues. 
This analysis will be carried out in a separate paper [To| . 

C. The quasifermion pole and quasifermion 
dispersion law 

1. How to define the quasifermion pole 

Generally speaking, the pole of the propagator or the 
point where the propagator inverse vanishes defines the 
corresponding particle and its dispersion law. In the 
case of the thermal quasi-particle, however, its mass term 
usually has finite, not small but in most occasion quite 
large imaginary part, namely, the pole position of ther- 
mal quasiparticle sits deeply inside the complex po plane. 

Because it is not very simple to study the structure of 
such a pole sitting deeply inside the complex po plane, 
we usually study such a pole by defining the condition 
so that the real part of propagator inverse vanishes as 
the onshell condition. We adopt this definition of onshell 
throughout this analysis, then the quasifermion pole is 
defined by the zero-point of the real part of chiral invari- 
ant fermion propagator inverse Re[D±(P) = D±(po,p)], 

Re[D±(p ,p)] = at (p = w±,p), (13) 

which determines the dispersion law of this pole, ui±(p). 

There is of course another definition of onshell and its 
corresponding pole. One of such definition is to use the 
peak position of spectral density as the pole position of 
the corresponding particle. With this definition we can 
also determine the dispersion law of this pole, u>±(j>). 
Though in most cases these two definitions give the same 
results, i.e., the dispersion law determined through the 



9 



onshell condition Ke[D±(po = lu±,p)} — agrees with 
the one determined by the peak position of the spectral 
density, in some cases two definitions give different re- 
sults. We have also discussed in the Sec. IIII Bl above, the 
possibility that the peak position of the spectral density 
in the region of intermediate coupling strength may not 
correctly represent the physical modes. On this problem 
we will discuss in Appendix |DJ Therefore, as mentioned 
above, we adopt Eq. ([13]) as the definition of onshell. 



2. Quasifermion dispersion law in the weakly coupled 
QCD/QED medium and the fermion thermal mass 

Now let us study the quasifermion dispersion law de- 
termined through the onshell condition, Eq. (|T3"|) . i.e., 
Re[D±(p — u±,p)] — 0. In Fig. [B] we give the 
quasifermion dispersion law uj = uj± (p) at small coupling 
and at moderately high temperature. It should be noted 
that, as can be seen in Fig. [5J in the region of weak cou- 
pling strength a < 0.01 the dispersion law lies on a uni- 
versal curve determined by the HTL calculations [ill (3 ■ 
Thus the result shows a good agreement with the HTL 
resummed effective perturbation calculation. 

The important point is that both the quasifermion en- 
ergy and the plasmino energy u>-(p) approach to 
the same fixed value rri*j, Eq. (I12[) . as p — > 0, namely, 
in Fig. [BJ the normalized energy w±(p) = uj±{p)/m^ ap- 
proach to 1 as p — > 0. This fact clearly shows that the 
quasifermion as well as the plasmino have a definite ther- 
mal mass m*f of O(gT) determined through the next- 
to-leading order calculation of HTL resummed effective 
perturbation theory 0, [HJ . We should also note that the 
collective plasmino mode exhibits a minimum at p ^ 
and vanishes rapidly on to the light cone as p gets large. 

It is also worth noticing that at weak coupling and 
moderately high temperature, dispersion law in the 
small-p region determined by the zero-point of D + agrees 
well and almost coincides with that determined by the 
peak of the spectral density p + . This fact can be under- 
stood by the one already noted in Sec IIII A"2l that in the 
weak coupling region the thermal mass of quasifermion 
determined by the peak position of spectral density al- 
most coincides with the one determined by the solution 
of the onshell condition Re[D + (p ,p = 0)] = 0, and that 
the thermal mass thus determined is m*j. The sharp 
peak structure of Im\B{po,p — 0)], or the narrow width 
structure of the quasifermion pole in the correspond- 
ing region may guarantee this fact. As the momentum 
p becomes large, however, there appears a discrepancy 
between them, especially in the dispersion law of the 
plsmino branch, see Appendix [D] 




FIG. 6: The normalized quasifermion dispersion law ui*(p) = 
u)(p)/m*f at T — 0.3 in the small coupling region. 



3. Vanishing of the thermal mass in the strongly coupled 
QCD/QED medium 

Next let us study how the result shown in Fig. [BJ 
changes as the coupling gets stronger, namely, in the 
region of intermediate to strong couplings. For this 
purpose, let us see carefully the fermion branch of the 
quasifermion dispersion law in the small momentum re- 
gion. (N.B. Temperatures and couplings we are studying 
belong to the chiral symmetric phase.) 

Figure [Jj shows the a-dependence of the normalized 
dispersion law at T = 0.3 as the coupling a becomes 
stronger, where the normalization scale is the next-to- 
leading order thermal mass m*^. In this Fig. [7] we can 
clearly see, though in the weak coupling region we get 
the solution in good agreement with the HTL resummed 
perturbation analyses, as the coupling becomes stronger 
from intermediate to strong coupling region the normal- 
ized thermal mass = 0) = uj + (p = 0)/my begins 
to decrease from 1 and finally tends to zero (a > 0.27 
in Fig. [7J. Namely, in the thermal QCD/QED medium, 
thermal mass of quasifermion begins to decrease as the 
strength of coupling gets stronger and finally disappears 
in the strong coupling region. This fact strongly suggests 
that in the recently produced strongly coupled QGP, 
thermal mass of quasifermion should vanish or at least 
become significantly lighter compared to the value in the 
ideal weakly coupled QGP. 

To see the above behavior of thermal mass more 
clearly, in Fig. [5] we show the normalized mass Lo+(p = 0) 
as a function of a. In the small coupling region (a < 0.1) 
and around the temperature T = 0.1 ~ 0.2, results of 
thermal mass agree well with those of the HTL resummed 
perturbation calculation. As the coupling gets stronger 
from intermediate to strong coupling regions, however, 
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FIG. 7: Thea-dependence of the normalized quasifermion dis- 



persion law u*(p) 
becomes stronger. 



uj(p)/m*f at T = 0.3, as the coupling 



the normalized thermal mass 0J+(p = 0) begins to de- 
crease from 1 and finally goes down to zero, i.e., the 
thermal mass vanishes. 
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FIG. 8: The a-dependence of the normalized thermal mass 
bJ+(j> = 0), see Text. 

Analogous behavior of thermal mass uj+(p = 0) ap- 
pears in the temperature-dependence. As can be seen in 
Fig. [51 with any coupling a, thermal mass decreases from 
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as the temperature becomes higher, and finally at 
extreme high temperature = 0) becomes zero, thus 
the thermal mass vanishes. Figure O shows another char- 
acteristic behavior as T — ¥ small. Almost at any coupling 
the thermal mass decreases and finally tends to vanish as 
temperature becomes lower. This behavior is consistent 
with the fact that at zero-temperature the thermal mass 
must vanish. The unexpected behavior is that, as the 
coupling becomes stronger, the thermal mass uj+(j> = 0) 
vanishes at low but non-zero finite temperature. 

Here it is to be noted that the ratio to+{p = 0)/T is not 
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FIG. 9: The a-dependence of the normalized thermal mass 
uj+(p = 0), see Text. 



necessarily a constant and the T-dependence at not-so- 
high T observed in our analysis is a consequence of the 
nonperturbative DSE analysis. It is because the addi- 
tional dimensionful parameter, such as the regularization 
(or the cutoff) scale or the renormalization scale comes 
into the theory through the regularization and/or the 
renormalization of massless thermal QCD/QED. In our 
case the cutoff scale A is introduced into the theory. The 
thermal mass, in fact, has a logarithmic T-dependence in 
the effective perturbation calculation, see, e.g., Rcbhan's 
lecture in Ref. [|[. 

The behavior of thermal mass is determined by the 
behavior of chiral invariant mass function Re[B(po,p)]. 
In Fig. [TU] we show, for the sake of convenience, the in- 
dependence of Re[D + (p ,p = 0)] =Re[po + B(p ,p = 0)] 
at T = 0.3. At small coupling Re[D + (p ,p = 0)] has a 
steep valley/peak structure in the small po region, but as 
the coupling becomes stronger this structure eventually 
disappears and Re[D + (p ,p = 0)] belongs to behavior 
almost as a straight line with a slope +1. 
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FIG. 10: The a-dependence of the real part of the inverse 
fermion propagator at p = 0, Re[D + (po,p = 0)]=Re[po + 
B(p ,P = 0), at T = O.3.] 

Thermal mass is given by the solution of Re[_D + (po,p = 
0)] = 0, i.e., po-coordinate of intersection point of the 
drawn curve of Re [D + (po ,p = ) ] and the po axis . At first 
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we can see with this figure that at small couplings there 
are three intersection points, the one with positive po, the 
one with negative po, and the one at the origin pa = 0, 
which correspond to the quasifcrmion, the plasmino and 
the massless (or, ultrasoft) modes 0,0], respectively. 

As the coupling becomes stronger (a > 0.27 at T = 
0.3), however, the number of the intersection points sud- 
denly reduces and there appears only one intersection 
point at po = 0, which may correspond to the massless 
pole in the fcrmion propagator. Thus we can understand 
the behavior in Fig. [FJ namely, in the weak coupling re- 
gion ui+ip = 0) = <jJ + (p = Q)/m*j is almost unity, and 
reduces to zero in the strong coupling region (a > 0.27 
at T — 0.3), showing that the fermion thermal mass van- 
ishes completely in the corresponding strong coupling re- 
gion. 



4- Disappearance of the plasmino mode in strongly coupled 
QCD/QED medium 

Finally we study what happens in the plasmino mode 
in the strongly coupled QCD/QED medium, by explic- 
itly examining the plasmino branch of the dispersion 
law. In the last Sec. IIII C 3[ where we see the thermal 
mass vanishes in the strong coupling region, we only 
studied the structure of the fermion branch of the dis- 
persion law, and of the propagator inverse at p = 0, 
Re[D + (p ,p = 0)] = 0. We can not exactly see what 
happens in the plasmino mode without explicitly study- 
ing the plasmino branch of the dispersion law. 

Figure [11] shows the a-dependence of the normalized 
dispersion law of the plasmino branch at T — 0.3 as the 
coupling a becomes stronger. (At weak couplings we 
already saw its structure in Fig.[BJ) 

Paying attention to the plasmino branch, we recognize 
that, as the coupling gets stronger, the valley structure 
of the plasmino dispersion law, or, the existence of the 
minimum in the plasmino dispersion law, observed in the 
weak coupling region, eventually disappears, and that the 
plasmino dispersion law gets to sharply dropping onto the 
light cone as the momentum p becomes large. 

At T = 0.3 in the small coupling region a < 0.01 the 
plasmino branch lies on the universal curve determined 
by the HTL calculation. Around a ~ 0.02 the dispersion 
law of the plasmino branch begins to change its structure: 
firstly the behavior as p — > large begins to show sud- 
den decrease onto the light cone, then secondly around 
a ~ 0.05 the valley structure of the plasmino dispersion 
law eventually disappears and the plasmino dispersion 
law monotonically drops sharply onto the light-cone, and 
finally in the region a > 0.27 (at T — 0.3) the plasmino 
branch totally disappears. 

If the coupling gets further stronger, the thermal mass 
begins to decrease and eventually disappears at a > 0.27, 
as noted in Sec. IIII C~3l before. The plasmino branch dis- 
appears at a > 0.27 also, which we can see in Fig. [TT1 
and the three modes, i.e., the fermion, the plasmino and 




p/m* f 

FIG. 11: The Q-dependence of the normalized dispersion law 
u>*L(p) = uj-(p)/m*f of the plasmino branch at T = 0.3. The 
curve shown is the dispersion law of the plasmino branch in 
the HTL calculation [lj 

the ultrasoft modes, finally merge and become a single 
massless mode that can be hardly detected as a real phys- 
ical mode in the strongly coupled QGP, as noted before 
because of its large decay width. 

D. Thermal mass of the quasifermion 

In the last Sec. IIII O we have disclosed unexpected be- 
havior of the thermal mass of quasifermion in the strong 
coupling QCD/QED, namely the fact that the ther- 
mal mass vanishes in the strongly coupled QCD/QED 
medium (or, the recently produced strongly coupled 
QGP). Also we have pointed out that at weak coupling 
and high temperature both the quasifermion and the 
plasmino have a common thermal mass m% , Eq. (|12[) . de- 
termined through the next-to-leading order calculation of 
HTL resummed effective perturbation theory [3j, |21| . 

In this section we examine how accurately the ther- 
mal mass ml Eq. (|12l) . can describe the thermal mass 
calculated in our analysis. Fig. [8j showing the coupling 
dependence of the thermal mass presented in the last sec- 
tion, covers a wide range of couplings and gives us only 
a rough image, thus is not suited to the present purpose. 

Here we present Fig. [TSJ the rescaled version of Fig. [51 
showing the thermal mass calculated in our analysis in 
the weak coupling region a < 0.1. Now we can see 
clearly that in the whole region of the temperature 
0.100 < T < 0.200 at weak couplings a < 0.1, the nor- 
malized thermal mass = 0) = u> + (p = fy/m** is 
almost unity, namely, the thermal mass calculated in our 
analysis, co+(p = 0), is well described by m*j. As the 
temperature becomes higher, discrepancy becomes evi- 
dent and larger; the normalized thermal mass lo+(j> = 0) 
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FIG. 12: The a-dependence of the normalized thermal mass 
aj+(p = 0), see Text. 

becomes to deviate from unity and gets smaller, namely, 
uj+(jp — 0) begins to decrease from rrA and becomes 
smaller. 

It should be noted, however, while at very small cou- 
plings a < 0.01 a common tendency can be recognized 
in Fig. [5] that the normalized thermal mass ui+(p = 0) 
approaches unity, except at extreme high temperatures. 

In studying the temperature dependence of the ther- 
mal mass, which is shown in Fig. |9l another fact can be 
recognized. The hrst thing that attracts our attention 
is that, except in the small coupling region, the normal- 
ized thermal mass = 0) shows a peak structure, 
namely that uj+(p = 0) decreases as both the temper- 
ature becomes higher and becomes lower. This fact, on 
the former case we already noticed in Sec. MI C 31 is unex- 
pected and not easy to be understood with the knowledge 
we have learned through the effective perturbation anal- 
yses. The behavior in the lower temperature region may 
indicate that the thermal mass shows a behavior pro- 
portional to T/log(l/T), while in the high temperature 
region the thermal mass shows a behavior proportional 
to Tlog(l/T). (N.B.: The temperature T varies in the 
range < T < 1.) 

As for the temperature dependence in the higher tem- 
perature region, only we can say definitely at present 
that the ratio ui+(p = 0)/T is not necessarily a constant 
and the T-dependence of the normalized thermal mass 
= 0) observed in our analysis is a consequence of 
the nonperturbative DSE analysis. 

At small couplings the normalized thermal mass 
uj+(p — 0) seems to approach unity as the temperature 
becomes lower, showing the well-known behavior of the 
thermal mass being proportional to the temperature T, 
and is easy to be understood. 



E. Existence of the third peak, or the ultrasoft 
mode 

The quasifermion and the plasmino modes are well un- 
derstood in the HTL resummed analyses, the latter be- 



ing the collective mode to appear in the thermal envi- 
ronment. What is the third peak? Is it nothing but 
convincing evidence of the existence of a massless or an 
ultrasoft mode? Is there any signature in our analysis? 

Existence of the massless or the ultrasoft fermionic 
mode has been suggested first in the one-loop calcula- 
tion [H when a fermion is coupled with a massive boson 
with mass m. The spectral function of the fermion gets to 
have a massless peak in addition to the normal fermion 
and the plasmino peaks. Recently a possible existence 
of collective fermionic excitation in the ultrasoft energy- 
momentum region p < g 2 T is investigated analytically 
through perturbative calculation [H, H4I . Both of these 
analyses are confined to the weak coupling regime, and 
nothing is known what happens in the sQGP we are in- 
terested in. In this sense first we will study the structure 
of the third mode, i.e., of the massless or the ultrasoft 
mode, in the weakly coupled QCD/QED medium, then 
will proceed to the intermediate and strong coupling re- 
gion to investigate how the ultrasoft mode behaves in 
such an environment [9(. 

Now let us study the structure of the third mode, i.e., 
of the massless or the ultrasoft mode, in the weakly cou- 
pled QCD/QED medium. First we give in Fig. flUTa) 
the structure of spectral density p + (po,p = 0) at weak 
coupling region (a = 0.001, T — 0.3). Two sharp peaks, 
representing the quasifermion and the plasmino poles are 
clearly seen, and the existence of a slight "peak" can also 
be recognized around po — 0. To see more clearly, in 
Fig. ITBT b') shown is rescaled version of Fig. [T37a). where 
we can clearly see the "peak" structure around po = 0. 
This third peak is nothing but a convincing evidence of 
the existence of a massless or an ultrasoft mode [J, |9|. 
This peak is indistinctively slight compared to the sharp 
quasifermion and plasmino peaks. 

Here we should take notice of the fact that the peak- 
height of the ultrasoft mode centered at po = is, roughly 
speaking, O(g) lower than the peak-height of the normal 
fermion or the plasmino peak centered at po = 
This rough result does not exactly agree with what Hi- 
daka, Satow and Kunihiro have shown in their works [22] 
concerning the residue of the ultrasoft fermion mode, and 
we will perform more detailed analysis on this problem 
in a separate paper. 



F. Decay width of the quasifermion, or the 
imaginary part of the chiral invariant mass function 

B 

Finally let us study the decay width of the 
quasifermion or the imaginary part of the chiral invariant 
mass function Im[B(po,p)] at p = 0. The decay width of 
the quasifermion is extensively studied through the HTL 
resummed effective perturbation calculation [241 ] ■ giving 
a gauge-invariant result of 0(g 2 T\og(l/g)). However, as 
is shown above, the quasiparticle exhibits an unexpected 
behavior, such as the vanishing of the thermal mass in the 
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FIG. 13: (a) Quasifermion spectral density p+(po,p = 0) at 
small coupling region (a = 0.001, T — 0.3). (b) Quasifermion 
spectral density p + (po,p = 0) enlarged around the origin. 



strongly coupled QCD/QED medium, completely differ- 
ent from that expected from the HTL resummed effective 
perturbation analyses. How does the decay width of the 
quasifermion exhibit its property in the corresponding 
strongly coupled QCD/QED medium? 
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FIG. 14: The decay width of the quasifermion at rest *y(p = 0) 
in the weakly coupled QGP. The dotted straight line Eq. (|15[) 
represents the result from the HTL resummed calculation, see 
Text. 



In Fig. H~4l and Fig. [15] we show the decay width of the 
quasifermion 7(p) at p = 0, in the weakly coupled and in 
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FIG. 15: The decay width of the quasifermion at rest f(p — 0) 
in the strongly coupled QGP. The horizontal straight line 
Eq. (|15|) represents the result from the HTL resummed cal- 
culation, see Text. 



the strongly coupled QGP, respectively, where 



lip) = 2 Im [- D +(^° = UJ +(P)>P)] 



—Rc[D + (p Q ,p)} . (14) 

°P° p =w+(p). 

In both figures, the fitting straight line represents 

lip = 0) = iaT^logi + c^) , (15) 
c ~ 3.38. 

In the weak coupling and high temperature QGP, the 
decay width j(p — 0), Eq. (fT5|) . agrees with the HTL 
resummed effective perturbation calculation [24| up to a 
numerical factor [25|], see Fig. Q3] 

Quite unexpectedly even in the strongly coupled QGP, 
the resulting decay width j(p = 0), Eq. (TL31) . shows the 
totally same behavior of 0{g 2 T log(l/s)) as in the weakly 
coupled QGP up to the numerical factor and the 0(g 2 T) 
correction term, Fig. 1151 

What happens in the intermediate coupling region? 
The results of the decay width at T — 0.150 is given 
in Fig. [16] From this figure we can understand how the 
decay width in the weak coupled QGP and the one in the 
strongly coupled QGP coincide. In the intermediate cou- 
pling region, the decay width of the quasifermion shows 
a "rich" structure. The decay width j(p = 0) diverges 
at the vanishing point of the thermal mass to+{p = 0), 
namely, the point where ui+{p = 0) first hits zero as the 
coupling changes, see Eq. (|14l) . 

This behavior is again not expected, because the 
quasifermion in the small coupling and high temperature 
QGP and the one in the strong coupling and high tem- 
perature QGP are totally different; in the former case the 
quasifermion has a thermal mass of O(gT) and the plas- 
mino branch exists in a fermion dispersion law, while in 
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FIG. 16: The a-dependence of the decay width of the 
quasifermion at rest j(p = 0) at T = 0.15. 

the latter case thermal mass of the quasifermion vanishes 
and the plasmino branch disappears. 

The temperature-dependence of the decay width is 
again described by Eq. (TToT) . namely, the decay width 
of the quasifermion is linearly proportional to the tem- 
perature, both in the weak and strong coupling QGP. 

This behavior can be clearly seen in Figs. [17] and [T8l 
and also in Figs. and [11)] The former fact indicates 
that in the strongly coupled QGP, recently produced at 
RHIC and LHC, the predicted massless or the ultrasoft 
pole is very hard to be detected as a real physical mode. 



0.9 

~ 0i8 a=0.01 O 

™ „ , a=0.02 

n a=0.03 o 

i „ „ a=0.04 O 

£ 0-6 a=0.05 

| 0.5 

| o.3 $ —«—&—•§ 6 'I t § — i! 

Q, 

1= 0.2 
0.1 

o I ' ' ' ' ' ' ' ' ' 

0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.2 

T 

FIG. 17: The T-dependence of the decay width of the 
quasifermion at rest ~f(p = 0) in the weakly coupled QGP. 
The dotted straight line Eq. (|15[) represents the result from 
the HTL resummed calculation, see Text. 



IV. SUMMARY AND DISCUSSION 

In this paper we carried out a nonperturbative anal- 
ysis of a thermal quasifermion in thermal QCD/QED 
by studying its self-energy function through the Dyson- 
Schwinger equation (DSE) with the HTL resummed im- 
proved ladder kernel. With the solution of the DSE 
we studied the properties of the thermal quasifermion 



FIG. 18: The T-dependence of the decay width of the 
quasifermion at rest j(p — 0) in the strongly coupled QGP. 
The horizontal straight line Eq. (|15[) represents the result from 
the HTL resummed calculation, see Text 

spectral density and its peak structure, the dispersion 
law of the physical modes corresponding to the poles of 
thermal quasifermion propagator. Through the study 
of quasifermion we elucidated the properties of ther- 
mal mass and the decay width of fermion and plasmino 
modes, and also paid attention to properties of the pos- 
sible third mode, both especially in the strongly coupled 
QCD/QED medium. 

What we have revealed in this paper is the drastic 
change of properties of the "quasifermion" depending on 
the strength of the interaction among constituents of the 
QCD/QED medium; 

i) In the weak coupling region, or in the weakly 
coupled QCD/QED medium: a < 0.02 or g < 0.5 
at T = 0.3. The onshell conditions through the 
peak structure of spectral density and from the 
zero point of the quasifermion inverse propaga- 
tor give the same structure and properties of the 
quasifermion. A rigid quasiparticle picture holds 
with the thermal mass mp Eq. (|12p . and a small 
imaginary part or the decay rate 7 ~ g 2 T 4og(l/ g) 
and the fermion acts as a basic degree of freedom 
of the medium. The thermal mass is nothing 
but the next-to-leading order result of the HTL re- 
summed effective perturbation calculations Q. A 
fermion and the plasmino mode appear. Thus the 
results in the weak coupling region well reproduce 
those of the HTL resummed effective perturbation 
calculations. 

The triple peak structure of the quasifermion spec- 
tral density is clearly observed, indicating the exis- 
tence of the fcrmionic ultrasoft third mode, which 
is absent from the HTL resummed effective pertur- 
bation analyses. 

ii) In the strong coupling region, or in the strongly 
coupled QCD/QED medium: a > 0.1 or g > 1 at 
T = 0.3. Both the spectral density and the inverse 
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fermion propagator tell the single massless peak 
structure with large imaginary part, or, the decay 
rate 7 ~ g 2 T\og(l/g). The quasiparticle picture a 
la Landau has been broken down in the strongly 
coupled QCD/QED medium. The thermal mass 
vanishes and there appears only the fermion mode, 
the plasmino mode disappears in the strongly cou- 
pled medium. 

hi) In the intermediate coupling region: 0.02 < a < 0.1 
or 0.5 < g < 1 at T = 0.3. In this region the spec- 
tral density and the inverse fermion propagator tell 
completely different structure. The spectral den- 
sity tells that there should be two particle modes 
with large decay rates, while the inverse fermion 
propagator tells that there should be three poles in 
the propagator, thus may exist three modes in this 
coupling region, just as in the case in the weakly 
coupled medium. We conclude that the indication 
of the inverse fermion propagator tells the truth, 
see the text. Anyway the intermediate coupling re- 
gion is the transitional region for the fermion in the 
medium to behave as a rigid quasiperticle, acting 
as a basic degree of freedom in the medium. 

Here we give several comments and discussion on the 
results of the present analysis. 

1) It is not a priori very clear which one really defines 
the onshellness of the physical particle and its dispersion 
law, the peak position of the quasifermion spectral den- 
sity p±(po,p) or the zero-point of the real part of the 
inverse fermion propagator ~Re[D±(po,p)] — 0, especially 
when the imaginary part is not very small. If we adopt 
the peak position of the quasifermion spectral density 
p±(po,p) as the onshell point of the particle, then the 
corresponding dispersion law exhibits a branch develop- 
ing into the space- like domain of space-time. There is also 
the problem of double peak structure of the quasifermion 
spectral density in the transitional intermediate coupling 
region, as noted above in iii). With these facts we adopt 
Ke[D±(po,p)] = as the onshell condition of the physical 
particle to study its dispersion law and various proper- 
ties, such as the thermal mass and the decay width, etc.. 

2) With the onshell condition Re[D±(po,p)] = we 
select the particle mode and study its dispersion law and 
the particle properties. The onshell condition at p = 0, 
Re[D±(p ,p = 0)] = 0, always has solution at po = 0, 
which may correspond to the ultrasoft mode. This cor- 
respondence is, however, not so simple. The structure of 
the imaginary part around the onshell point of the prop- 
agator plays an important role to make this correspon- 
dence exact. This relationship was pointed out by Ki- 
tazawa et al. jj| ; if in the medium the bosonic mode with 
nonzero mass (with small decay rate) couples with the 
fermion, then the quasifermion spectral density shows a 
triple peak structure corresponding to the ultrasoft third 
mode together with the fermion and the plasmino modes. 
The appearance of the peak at po = corresponding to 
the ultrasoft third mode is guaranteed with the vanish- 



ment of the imaginary part, Im[D±(pQ,p = 0)] = 0, at 
Po = 0, which happens because of the coupling of the 
fermion with the massive bosonic mode. Such a mecha- 
nism may not work in the QED medium since no massive 
bosonic excitations in the QED medium are expected. 

In the present DSE analysis, the imaginary part of 
the fermion propagator does not vanish at po = 0, 
Im[D±(po,p = 0)] 7^ 0, but rather shows a peak struc- 
ture at po = 0. This peak structure actually suppreses 
the peak-height of the spectral density, as noted in the 
text, Sec. lIIIEl In this sense the appearance of the ultra- 
soft third peak with very low peak-height in the present 
analysis may have different origin from that of Kitazawa 
[j| and from Satow et al. 



et 



22] 



roblem will 



Thisp: 

be discussed further in a separate paper [10(. 

3) We have noted that the thermal mass of 
quasifermion decreases as the coupling gets stronger, and 
finally vanishes in the strong coupling region. This fact 
indicates that the thermal mass of quasiquark vanishes 
and behaves as massless fermion with large decay rate in 
the recently discovered strongly coupled QGP. It is not 
so simple, however, whether such a particle can be ex- 
perimentally observed as a massless quasifermion or not. 

4) As noted above, the decay rate j(p) of the 
quasifermion in the QCD/QED medium shows a typi- 
cal g 2 T\og(l/g) behavior both in the weakly and the 
strongly coupled medium. In the transitional intermedi- 
ate coupling environment, however, the decay rate shows 
a rich structure and j(p = 0) even diverges at the vanish- 
ing point of the thermal mass. It is quite exciting if we 
can find some methods to be able to verify experimen- 
tally the unexpected behavior of the thermal mass and 
the decay rate. 



Appendix A: Approximations to get the HTL 
resummed improved ladder DSE 

In the present analysis, we solve the DS equation for 
the retarded fermion self-energy function E#, with the 
HTL resummed gauge boson propagator Eq. ©, by 
adopting further the following two approximations to get 
Eq.|H i) the point- vertex approximation, and ii) the mod- 
ified instantaneous exchange approximation to get the fi- 
nal DSEs we solve, on which we give brief explanations 
below. 

i) Point- vertex approximation to get Eq. (j4]) 
As for the vertex function we adopt the point- 
vertex approximation, namely we simply set *F Al = 7^ 
disregarding the HTL corrections to *F'\ Thus we in- 
vestigate the ladder (point- vertex) DS equation with the 
HTL resummed gauge boson propagator. 

There are two reasons; Firstly, without the point- 
vertex approximation the numerical calculation we 
should carry out becomes so complicated that we can not 
manage with the power of the computer we use, because 
the HTL resummed contribution to the vertex function 
is the non-local interaction term, and also because it be- 
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haves singular in numerical calculations. Secondly, in the 
DS equation with the HTL resummed vertex function, it 
is difficult to resolve the problem of double counting of di- 
agrams [26j], especially at the level of numerical analyses. 
Being free from this problem in the numerical analysis 
is the main reason why we make use of the point-vertex 
approximation. 

ii) Modified Instantaneous Exchange (MIE) approxi- 



mation to get the final DSEs to solve 

The next approximation we make use of is the modified 
instantaneous exchange (IE) approximation (i.e., set the 
energy component of the gauge boson to be zero) to the 
gauge boson propagator *G fJ/U . The retarded (R = RA) 
and correlation (C = RR) components of the HTL re- 
summed gauge boson propagator *G^ V are given by [l6j 



*G%(K) = *G% A {-K, K) 



*G%{K) = *G% R {-K,K) 



Ji^v , 

*H$(K) - K 2 - iek 

(l + 2n B (k ))[*G R »(K)- 



*n£ (K) - k 2 

'G%(K)] t 



iekn 



i 

K 2 + iek Q 



D^, (Ala) 
(Alb) 



with 



"11% and 



TI^ being the HTL contributions to the 



transverse and logitudinal modes of the retarded gauge 
boson self-energy, respectively [15j. The parameter £ is 
the gauge-fixing parameter (£ = in the Landau gauge). 

In the above, A^, B^ v and are the projection 
tensors given by [l6j 



9" 



K^K h 
K 2 " 



(A2a) 
(A2b) 
(A2c) 



where K = (fc, kok), k = vk 2 and k = k/fc denotes the 
unit three vector along k. 

The modified IE approximation we make use of con- 
sists of taking the IE limit in the HTL resummed logi- 
tudinal (electric) gauge boson propagator, *G^", that is 
proportional to B^ v , while keeping the exact HTL re- 
summed form for the travsverse (magnetic) gauge boson 
propagator, *Gt^ , that is proportional to A^ v , and also 
for the massless gauge term in proportion to . The 
reason why we do not take the IE limit to the transverse 
mode is that the IE approximation reduces the transverse 
mode to the pure massless propagation, thus makes the 
important thermal effect, i.e., the dynamical screening of 
transverse propagation disappear. 

With the above two approximations, we obtain the fi- 
nal HTL resummed improved ladder DS equations for 
ther invarinat scalar functions A, B and C, to solve. 



Appendix B: Cutoff dependence 

In this appendix we explain the cutoff dependence of 
the present analysis. 

As explained in Sec. Ill A[ in solving the DSEs, Eq. ((U), 
we are forced to introduce a momentum cutoff in the in- 



tegration over the four-momentum J d 4 K, K = (fco,k) 
is the fermion four-momentum. The cutoff method we 
make use of is as follows (A denotes an arbitrary cut- 
off parameter and plays a role to scale any dimensionful 
quantity, e.g., T = 0.3 means T = 0.3A); 



three momentum k 
energy fc 



k = |k| < A 
1*3 1 < A 



In the present analysis we make the ratio r = Aq /A vary 
r = 1 ~ 5, and fix it so as to get a stable solution to the 
fermion spectral density. 

In Fig. [T2] we show how the spectral density changes 
as a function of kg as we vary the ratio in the range 
r = 1 ~ 5. We can easily recognize that at T = 0.4 and 
a = 0.1 we can get a stable solution if we choose r > 2. 

Situation is almost the same but slightly differs at dif- 
ferent T and a, see Fig. HOI at T = 0.4, a = 1.0 and 
compare with Fig. Q2] at T = 0.4, a — 0.1. As the cou- 
pling becomes stronger it is safe to choose larger values 
of r. 

The stability of the solution can be checked by the 
saturation of the sum rules, Eqs. (|10aj) and (|10b[) . As 
already noted, the sum rule Eq. (I10c[) heavily relied on 
the HTL calculation, thus we do not use this sum rule. 
Result is given in Table HI again showing the stability of 
the solution when we choose r > 2 (or more safely r > 3). 

With the above results we choose, in most cases ex- 
cept at very strong couplings, an appropriate value of 
the ratio in the range r > 2, depending on the region of 
temperatures and/or couplings we study. The extreme 
high temperature may cause another problem, namely, 
the problem of simulation artifact, therefore we restrict 
the temperature to the region T < 0.6 and do not per- 
form our analysis in an extreme high temperature region 
T > 0.8. 
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FIG. 19: The cutoff dependence of the quasifermion spectral 
density at p = 0, p+(po,p ~ 0), at T = 0.4 and a — 0.1. 
r = Aq/A. 



FIG. 20: The cutoff dependence of the quasifermion spectral 
density at p = 0, p+(po,p = 0), at T = 0.4 and a — 1.0. 
r = Ao/A. 



TABLE I: The cutoff dependence of the saturation of the sum rules, Eqs. (10a) and (10b): r = Aq/A. 





r = 1 


r = 2 


r = 3 


r = 4 


r = 5 


T = 0.2 
a = 0.01 


Eq. ifTOa)) 


p = 


0.949 


0.994 


0.994 


0.996 


1.003 


p = 0.02 


0.986 


0.985 


0.983 


0.982 


0.990 


Eq. (fT0b| 


p = 0.02 


0.0197 


0.0207 


0.0207 


0.0209 


0.0202 


p = 0.04 


0.0392 


0.0404 


0.0404 


0.0407 


0.0404 


p = 0.1 


0.0990 


0.1009 


0.1009 


0.1004 


0.1013 


T = 0.4 
a = 0.1 


Eq. (fTOa)) 


p = 


0.927 


1.006 


1.006 


1.006 


1.006 


p = 0.02 


0.929 


1.006 


1.006 


1.006 


1.006 


Eq. (fTObl 


p = 0.02 


0.0180 


0.0206 


0.0207 


0.0207 


0.0207 


p = 0.04 


0.0354 


0.0406 


0.0408 


0.0408 


0.0408 


p = 0.1 


0.0948 


0.1011 


0.1016 


0.1016 


0.1016 


T = 0.4 
a = 1.0 


Eq. (fTOa)) 


p = 


0.759 


0.974 


1.018 


1.021 


1.021 


p = 0.02 


0.761 


0.975 


1.017 


1.021 


1.021 


Eq. (fTObl 


p = 0.02 


0.0108 


0.0171 


0.0199 


0.0202 


0.0203 


p = 0.04 


0.0218 


0.0350 


0.0407 


0.0413 


0.0413 


p = 0.1 


0.0546 


0.0879 


0.1024 


0.1037 


0.1038 



Appendix C: Phase boundary in the Landau gauge 

In order to study the phase transition and to determine 
the phase boundary of thermal QCD/QED, we should 
solve the DSE for the retarded fermion self-energy func- 
tion Sfl, Eq. @. For the present purpose, however, we 
must study the that has a c-number scalar mass func- 



tion C(P), 

E fl (P) = (1 - A(P))piY - S(F) 7 ° - C{P). (CI) 

The DSE in the Landau gauge to determine the three 
scalar invariants A(P), B(P) and C(P) becomes coupled 
integral equations as follows; 
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P 2 [1-A(P)} = g 



d 4 K 



(2tt)4 
Pp + 



{1 + 2n B ( PQ - fc )}Im[ *G%(P - K)} x 



{K a P p + K p P a - po{K a g pQ + K p g aQ ) - fc (P -.g P o + P p g a a) + pkzg ap 

A(K) 



-2p g a0 g p0 }- 



[ko + B{K) + ie} 2 - A(K) 2 k 2 - C(K) 2 
k + B(K) 



{p<t9 p q + P P g, 



<t() 



+ {l-2n F (k )} x 



[Ad + B[K) + ie} 2 - A(K) 2 k 2 - C(K) 2 . 
*G%(P - K)lm[{K a P p + K p P a - p {K a g p0 + K p g^) - ko{P a g p o + P p g M ) 

A(K) 



+pkzg ap + 2p k g <TO g p o} 



B(P) = g 



+{Pa9 P o + P P 9aO - 2pog a0 g p0 } 
d A K 



[ko + B(K) + ie] 2 - A(K) 2 k 2 - C(K) 2 
k + B(K) 



[k^ + B{K) + ie] 2 - A{K) 2 k 2 - C{K) 2 
{1 + 2n B {p^ - fc )}Im[ *G%(P - K)] x 

A(K) 



(2^) 4 

{K a g P o + K p g a0 - 2k g a0 g p0 } 



+{2g P a2g a0 - g ap \ 



*G p B a (P~ K)lm 



[k + B{K) + ie] 2 - A(K) 2 k 2 - C(K) 2 
k + B(K) 



[ko + B(K) + ie\ 2 - A(K) 2 k 2 - C(K) 2 
A(K) 



-2Ao5ao5po} 



[ko + B(K) + ie} 2 - A(K) 2 k 2 - C{K) 2 
k Q + B(K) 



+ {1 - 2n F (k )} x 
{K a g P o + K p g a o 



[ko + B(K) + ie] 2 - A{K) 2 k 2 - C(K)< 



{2g p o2g a0 - 9a P } 



C(P) = 9 2 



d 4 K 

(27T)" 



-9po 



{1 + 2n B {po - fc )}Im[ *G%(P - K)\ x 
C{K) 



[ko + B{K) + ie} 2 - A(K) 2 k 2 - C(K) 2 



+ {l-2n F (fc )} x 



f G p °(P-K)Im 



C(K) 



[k + B(K) + ie} 2 - A(K) 2 k 2 - C(K) 2 



(C2a) 



(C2b) 



(C2c) 



Above DSEs, Eq. (|C2j) . may have several solutions, and 
we choose the "true" solution by evaluating the effective 
potential V[Sji] for the fermion propagator function Sr, 
then finding the lowest energy solution. The effective 
potential V[S#] we evaluate, is given in Sec. Ill Bl Eq. (|5|). 

Now we present Fig. showing the phase boundary 
curve in (T, 1 / a)-plane in Landau gauge, that separates 
the chiral symmetric phase from the broken one. This 
critical curve shows that the critical coupling inverse l/a c 
is a monotonically decreasing function of the temperature 
T slightly concave upwords, and displays two character- 
istic behaviors: (1) as T becomes lower, the critical cou- 
pling inverse \/ot c becomes larger and seems to increase 
from below to the zero temperature value [l3j], and 
(2) the critical temperature T c increases as the coupling 
inverse becomes smaller (coupling become stronger) , with 
possible saturation behavior approaching To ~ 0.21 from 
below in the strong coupling limit. 

It is worth noting that the critical coupling at zero 



temperature is af S {T = 0) = ^-n ~ 1.152 [Hj], and our 
result predicts slightly larger critical coupling a c (T — 
0) = 1.32 ±0.14. The errors are given in the 3a accuracy 
level of the least x 2 fit. Phase transition occurs only in 
the region a > a c (T = 0) ~ 1.32 and T < T ~ 0.21, 
i.e., chiral symmetry broken phase is restricted to the 
region of the (T, l/a)-plane lower than the critical curve 
in Fig. [2TJ Therefore it is obvious that the region of the 
coupling and temperature where we study the property 
of the quasifermion, is well inside the chiral symmetric 
phase. 



How does the property of the quasifermion change in- 
side the chiral symmetry broken phase? This is an inter- 
esting question. Does the quasifermion mode still exist 
in the broken phase? These questions will be discussed 
in a separate paper. 
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FIG. 21: The phase boundary curve in the (T, 1/a) plane 
in the Landau gauge, separating the chiral symmetric phase 
from the broken one. The error-bar assigned to the best fit 
curve at T — denotes the error in the 3a accuracy level. 

Appendix D: Dispersion law us p (j>) determined 
through the peak position of the spectral density p+ 

Throughout this paper we determined the dispersion 
law of the thermal quasifermion with the onshell con- 
dition Re[D + (po,p — 0)] = 0. As was explained in 
Sec. IIII C[ generally speaking, the pole of the propaga- 
tor or the point where the propagator inverse vanishes 
defines the corresponding particle and its dispersion law, 
and we can use another definition of onshellness. One of 
such definition is to use the peak position of the spectral 
density as the pole position of the corresponding particle, 
with which we can also determine the dispersion law of 
this particle. 

These two definitions of onshellness almost agree with 
each other when the imaginary part of the mass term is 
small. In fact, in the weak coupling region at high tem- 
perature, the fermion branch of the dispersion law de- 
termined through the peak position of the spectral den- 
sity p+ almost completely coincides with the dispersion 
law, Fig. [BJ determined through the onshell condition 
Re[D + (p 0) P = 0)]=0. 

There are mainly two reasons why we adopt the on- 
shell condition Re[D + (po,p = 0)] = rather than that 
given by the peak position of the spectral density in the 
present analysis. The first reason is already explained 
in Sec. IIII A 21 It is pointed out there that at the inter- 
mediate coupling strength the spectral density exhibits a 
typical double peak structure, indicating the existence of 
two poles in the quasifermion propagator. This is, how- 
ever, not the case. There are actually three poles in the 
propagator in the corresponding coupling region. The 
third peak representing the ultrasoft third pole is com- 
pletely hidden under the big tails of the broad two peaks, 
thus cannot be observed. The position of the two peaks 
does not exactly represent the true position of the pole 
neither. These facts indicate that information obtained 
through the analysis of the spectral density itself is not 



complete but even misunderstanding. 

The second reason why did not adopt simply the peak 
position of the spectral density as the pole position of 
the corresponding particle, is the appearance of the plas- 
mino branch continuing to exist in the space-like do- 
main, namely the existence of the space-like plasmino 
solution. In Fig. (2D we show the dispersion law i>J p (p) de- 
termined through the peak position of the spectral den- 
sity which exactly corresponding to Fig. [SJ showing 




p/m* f 

FIG. 22: The normalized quasifermion dispersion law oj^/mj 
determined through the peak position of the spectral density 
at T = 0.3 in the small coupling region, see Text. 

the dispersion law determined through the onshell condi- 
tion Ke[D + (po,p — 0)] =0. Though the fermion branch 
almost completely agrees with each other, the plasmino 
branch exhibits a typical difference. In Fig. [51 the plas- 
mino branch exhibits a minimum at p and vanishes 
rapidly on to the light cone as p gets large. In Fig. |2"21 
the plasmino branch also exhibits a minimum at p =/= 0, 
and approaches rapidly to the light cone, then crosses 
the light cone and continues to exists in the space-like 
domain of the world sheet. 

With these two reasons, in the present analysis we do 
not adopt to define the onshell condition through the 
peak position of the spectral density. 
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